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Abstract 


A  standard  simplified  model  of  the  problem  of  determining  the  free 
electron  distribution  in  an  ionized  gas  by  radio  sounding  experiments  is 
discussed,  and  a  relation  between  the  elements  of  the  corresponding  scattering 
matrix  is  derived.  Under  the  conditions  assumed,  a  closed  fonn  expression 
for  the  charge  distribution  in  terras  of  a  reflection  coefficient  which  is  a 
rational  function  of  frequency  is  obtained. 
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1,   Int  rod-get  ion 

This  paper  will  deal  vd.th  certain  mathematical  aspects  of  the  radio 
sounding  method  for  investigating  the  negative  free  charge  density  of  an 
ionized  gas  such  as  the  ionosphere. 

A  typical  radio  sounding  experiment  is  conducted  by  sending  a  pulse 
modulated  electromagnetic  wave  into  the  gas.  The  time  delay  before  the  reflected 
pulse  returns  is  recorded  for  every  value  of  the  carrier  frequency  over  a  con- 
tinuous band  of  frequencies,  providing  data  which  then  become  available  for 
studying  the  effective  index  of  refraction  of  the  ionized  medium.  '-  -' 

The  index  of  refraction  presented  to  a  propagating  electromagnetic  wave 
by  an  ionized  gas  is  determined  for  all  practical  purposes  by  the  free  electrons 
in  the  gas.  The  square  of  the  index  of  refraction  n  is  given  approximately  by 
the  formula  L  -l 

(1.1)  n2  =  1  -  N«' 


mw  e_ 

0 

where  the  quantities  appearing  in  (1.1 )  are  defined  as  follows: 

N  is  the  density  of  the  free  electrons 
e  is  the  charge  of  a  single  electron 
m  is  the  mass  of  a  single  electron 
CO  is  the  angular  frequency  of  the  propagating  wave 

e  is  the  dielectric  constant  of  free  space. 

2 
The  expression  (l.l)  for  n  is  valid  only  to  the  extent  that  one  can  neglect 

losses  due  to  collisions  of  the  electrons.  If  the  mean  free  path  of  the  free 

electrons  is  very  long  COTipared  to  the  wavelength  of  the  propagating  wave,  the 

collisional  frequency  and  thus  the  losses  due  to  absorption  of  electromagnetic 

energy  by  the  gas  will  be  small. 
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It  will  be  assumed  here  that  the  electron  density  in  (1.1  )  is  a  function 
of  one  variable  x,  the  height,  i.e.,  that  the  electron  distribution  is  uniform 
ill  any  plane  orthogonal  to  the  direction  of  the  x  axis  in  a  Cartesian 
coordinate  system  imbedded  in  the  gas.  Further,  the  coordinate  system  will 
be  so  chosen  that  the  ionized  medium  occurs  only  in  the  region  where  x  is 
positive,  so  that 

(1.2)  N  =  0   when   x  <  0  . 

Given  all  the  assumptions  presented  so  far,  a  steady  state  periodic  electric 
field  u(x,k)  which  propagates  in  the  ionised  medium  must  s?.tisfy  a  differential 
equation  of  the  form 

(1.3)  ^  di^'^^     ^  r^  "  '^(^^p^-*'*^)  "   0  »    ("^M   =  ^  '"'h^"  X  <  0)  , 

where  k  is  the  wave  number  and  V(x)  is  proportional  to  the  electron  density  at 
X  but  is  independent  of  k. 

The  time  delay  data  of  a  radio  sounding  experiment  can  he  used  to  find 
the  phase  of  the  complex  reflection  coefficient  at  x  =  0  for  every  frequency 
within  the  experimental  band  LJ»L  J.  xhe  phase  of  the  reflection  coefficient 
at  other  frequencies  and  its  modulus  are  not  given  by  the  experimental  data, 
although  in  order  to  obtain  the  electron  distribution  it  is  necessary  to  know 
the  complete  reflection  coefficient  at  every  value  of  the  f requency L J » L  J ^ 
Therefore,  it  will  be  assumed  that  the  reflection  coefficient  is  given  completely 
for  all  frequencies  by  .some  kind  of  extrapolation  from  the  experimental  data. 
Since  these  data  are  not  sufficient  to  provide  all  the  necessary  information, 
it  should  be  expected  that  only  certain  features  of  the  resulting  electron 
density  function  will  actually  be  determined  correctly.  The  question  of  what 
features  of  the  electron  density  are  determined  correctly  by  means  of  the 
extrapolated  reflection  coefficient  is  an  important  but  difficult  one;  no 
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attempt  will  be  made  to  settle  it  here.  It  must  be  observed,  however,  that 
the  results  derived  from  the  extrapolated  reflection  coefficient  are  alwaj's 
at  least  as  reliable  as  any  result  obtained  from  the  given  data  alone.* 

In  the  sections  which  follow,  a  method  \-nll   be  given  for  reconstructing 
the  electron  distribution  from  a  knowledge  of  the  complete  reflection  coef- 
ficient at  all  frequencies.  It  will  be  seen,  in  fact,  that  the  result  can  be 
stated  explicitly  in  closed  form  when  the  reflection  coefficient  is  a  rational 
function  of  the  frequency.    The  procedure  for  extrapolating  to  the  complete 
reflection  coefficient,  starting  from  time  delay  data,  will  be  described  in  a 
future  article  which  will  deal  with  the  application  of  results  presented  here. 

Section  2  of  this  paper  can  be  considered  an  extension  of  the  introduction, 
for  nothing  which  is  discussed  there  is  new.  Some  of  the  general  well-known 
properties  of  the  scattering  matrix  associated  with  the  differentia]  equation 
Cl.3)  are  described  in  Section  2  for  convenient  reference  later  on.  Specific 
reference  to  the  literature  in  that  discussion  will  not  be  made  because  most 
of  the  material  comes  from  a  commonly-held  body  of  knowledge  difficult  to 
attribute  to  any  individual. 

In  Section  3  a  derivation  will  be  given  of  a  method  for  reconstructing  V(x) 
in  the  differential  equation  (1.3)  from  the  complete  reflection  coefficient 
at  x  =  0.  Although  the  method  has  already  been  presented  in  other  places  M>  K 
the  derivation  given  here  is  more  physically  motivated  and  slightly  more  general 
than  the  one  previously  used. 


While  the  standard  method  for  reconstructing  the  electron  distribution 
apparently  makes  use  of  nothing  but  the  available  time  delay  data  (See  ril,r6l) 
it  should  be  pointed  out  that  some  unjustified  hidden  assumptions  are  contained' 
in  the  procedure  (See  [2j,[7j.). 

It  was  found  by  Barghann  ■- -J  that  one  can  reconstruct  the  potential  for  the 
radial  Schrbdinger  equation  of  quantum  mechanics  from  a  knowledge  of  the 
asymptotic  phase  shift  of  the  wave  function  under  similar  circumstances. 


Section  U  is  a  discussion  of  the  scattering  matrix  elements  when  the 
reflection  coefficient  is  a  rational  function  of  the  wave  number  k. 

Section  5  is  devoted  to  oroving  a  theorem  connecting  the  asymptotic  form 
of  the  wave  u(x, k)  at  large  distances  with  the  reflection  coefficient  at  x  =  0 
under  the  assumption  that  the  reflection  coefficient  is  a  rational  function  of 
the  wave  number.  At  the  same  time  an  exolicit  closed  form  expression  for  V(x) 
in  terms  of  elementary  functions  is  derived.   It  will  be  seen  also  that  uCx^k'', 
itself,  can  be  given  explicitly  in  terms  of  elementary  functions. 

2.   The  Scattering  Matrix 

If  V(x)  in  the  differential  equation  (1.3)  satisfies  appropriate  conditions, 
any  solution  u(x,k)  has  the  asymptotic  behavior 

lim    Uy^   +  B^e"^^''  -u(x,k)]  =  0 
x  — ■>+oo  1  I 


(2.1) 


X  — >-00 


lim    jA_e^^"  +  B_e~^^^  -u(x,k)l  =  0  , 


where  the  coefficients  A+  ,  B±  are  independent  of  x  but  may  depend  on  k, 

Accordi.ng  to  (2.1)  the  asymptotic  forms  of  u(x, k)  at  large  positive  and 
negative  values  of  x  can  be  regarded  as  simple  waves  moving  to  the  left  and 
right  in  free  space  where  the  function  V(x)  vanishes.  On  the  basis  of  the 

convention  that  the  suppressed  periodic  time  factor  in  the  waves  is  e     a 

-ikx 
left  moving  wave  can  be  defined  as  one  a3:\'mptotically  proportional  to  e     and 

ikx 
a  right  moving  wave  as  one  asymptotically  proportional  to  e    .In  addition, 

an  incoming  wave  can  be  defined  as  a  right  moving  wave  at  -oo  or  a  left  moving 

wave  at  +oo ,  while  an  outgoing  wave  can  be  defined  as  a  right  moving  wave  at 

+C0  or  a  left  moving  wave  at  -oo  . 


The  set  o^  all  Incoming  waves  can  be  regarded  as  a  two  dimensional  vector 
space  of  column  vectors  such  that  the  first  component  of  each  vector  is  the 
amplitude  of  the  right  moving  part  and  the  second  component  is  the  amplitude 
of  the  left  moving  part.  Similarly,  the  outgoing  waves  can  be  represented  by 
a  two  dimensional  column  vector  space  for  which  the  ■first  component  of  each 
vector  is  the  amplitude  of  the  right  mo^/ing  part  and  the  second  component  is 
the  am.plitude  of  the  left  moving  part. 

Corresponding  to  every  differential  equation  of  the  form  (1,3)  whose 
solutions  have  the  as:'/mptotic  behavior  (2.1)  there  exists  a  matrix,  called  the 
scattering  matrJx,  which  transforms  each  vector  in  the  incoming  wave  vector 
space  into  a  vector  cf  the  outgoing  wave  vector  space.  If  the  incoming  wave 
vector  contains  the  amplitudes  of  the  incoming  part  of  u(x, k)  in  the  asymptotic 
limits  (2.1),  the  corresponding  outgoing  wave  vector  obtained  by  means  of  the 
scattering  matrix  will  contain  the  amplitudes  of  the  outgoing  part.  The 
scattering  matrix  has  the  form: 

/t(k)  ,  p(k)\ 
(2.2)  S  =/  j  J 

\r(k)  ,  xi^)J 

the  elements  of  S  are  called: 

the  transmission  coefficient  on  the  right:  t(k) 
the  transmission  coefficient  on  the  left  :  f  (k) 
the  reflection  coefficient  on  the  right  :  p(k) 
the  reflection  coefficient  on  the  left   :  r(k) 

If,  for  example, 

11m    ,u(x,k)  -  {k_e^^   +  B_e'^^)L  •=  0  , 
X  — >-oo 

and 

lim    |u(x,k)  -  (A^e^'^^  +  B^e"^^)^  =  0  , 

X  — >+Q0 
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thsn 


=  S 


or  eqiilvalently. 


(2.3) 


A^  -  A_t(k)   +  B^p(k) 
B_  •»  A_r(k)   +  B^xi^)   . 


Several  facts  aboiit  the  scattering  matrix  are  well  knovm.     If  V(x!*  in 
(1.3)  is  sufficiently  veil -behaved,   the  elements  of  S  are  analytic  function? 
of  the  -t'^ve  numher  k  in  the  complex  plane,   and  these  func+ions  are  regular  on 
the  real  axis.     Furthermore,   under  the  condition  on  V(x)  in   (1.3)  r(k)  and 
t(k)  are  reg\ilar  in  the  upper  half  plane  and  on  the  real  axis,   and  in  the  upper 
half  plane  they  obey  the  asymptotic  conditions   for  large  k: 

(2.b)  r(k)  =  O(^),    t(k)  =  1  +  O(^)    . 

The  following  relations  hold  in  general  for  real  k: 

(2.5)  r(-k)-?TF5-,  t(-k)-rn^,  p(-k)=^Tn,  ri'kh  KkT  , 

where  the  bar  over  a  number  indicates  the  complex  conjugate.   It  is  also 
knot-m  that  S  is  a  unitary  matrix  when  k  is  real? 

(2.6)  S*(k)  =  S"^(k)  ,  k  real. 
In  addition,  the  reciprocity  law  holds? 

(2.7)  t(k)  ^-n'k) 
for  all  k. 

Relations  (2.5),  (2.6)  and  (2.7)  imply  for  real  k; 

|r(k)|^+|t(k)|^  =  r(k)r(-k)+t(k)t(-k)  =■  1  (energy  conser^/ation) 

(2.8) 

p(k)rnrr+rnrjl(k)  -  p(k)t(-k)+r(-k)t(k)  -  O  (phase  law)  . 
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3,       The  Inverse  Problem 

Associated  with  the  differential  equation  of  (1.3)  was  the  condition 
(3.1)  V(x)  =0,  X  <  0  . 

Because  of  (3.1)  the  asymptotic  limit  of  u(x,k:)  given  by  (2.1)  for  x  ^>  -oo 

is  identical  with  u(x,k)  itself  for  negative  values  of  x.   If  in  (2.3) 

A_  =1,  B^  »  0  , 
then 

u(x,k)  -  e^^-  +  r(k)e"^^  ,  x<  0 
(3.2) 

and    lim    Ju(x,k)  -  t(k)e^^l  =  0  . 
X  — >  00  1  J 

The  asymptotic  conditions  (3.2)  mean  that  a  wave  of  unit  amplitude  is 
incident  from  the  left,  producing  a  reflected  wave  of  smplitude  r(k)  and  a 
transmitted  wave  of  amplitude  t(k).  Assuming  that  r(k)  is  given  for  all 
positive  real  values  of  k,  the  central  problem  of  this  paper  is  to  find  the 
function  V(x)  in  the  differential  equation  (1.3). 

In  accordance  with  the  physinal  problem  described  in  Section  1  it  is 
assumed  that  no  proper  eigenf unctions  of  (1.3)  e^xist.   Because  of  this  fact 
and  (3.1),  r(k)  is  regular  in  the  upper  half  plane '-^^-' . 

The  Fourier  transform  U(x,t)  of  u(x, k)  is  a  transient  wave  solution  of 
the  time  dependent  differential  equation 

(3.3)  B^VU^  _  £lT^  _  ^^^)  ^(^^^j  ^  0  ^ 

where  the  wave  velocity  has  been  taken  to  be  one.  If  u(x, k)  satisfies  the 

asymptotic  conditions  (3*2),  the  Fourier  transform  of  u(x,k)  exists  only  as 

as  improper  function;  i.e.,  it  contains  a  Dirac  delta  function  in  the  region 

x  <  0.  However,  if  the  incident  continuous  wave  of  unit  amplitude  is  replaced 

by  one  having  an  integrable  amplitude  ir  the  interval  (-od,+oo),  the  corresponding 


Since  the  charge  density  is  all  negative  V(x)  >  0  for  all  x. 
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transient  solution  of    (3.3)  is  a  proper  function.     The  transient  solution 
of   (3.3)   corresponding  to  an  incident  steady-state  continuous  v.'ave  of  unit 
ajriplitude  may  therefore  be  regarded  as  a  limiting  case  of  well-behaved 
transients  -  -  u&eful  because  of  its  formal  sinrlicity. 
The  reflected  transient  due  to  the  incident  wave 

inc  \        '   >  > 

is 

(3.)4)  R(xH)   =  jl     f         r(k)Q-^*'^^^*^  dk  . 

•^ -00 

Tn  order  to  obtain  (3»b)  vThen  r(k)  is  known  for  positive  real  k  it  is  neces- 
sary to  extend  r(k)  to  the  negative  real  axis  by  using  (2.5).  At  the  same 
time  it  can  be  seen  that  R(x+t)  is  real  for  all  real  x  and  t.  The  total  solu- 
tion U(x,t)  of  (3.3)  in  the  region  x  <  0  is  then 

(3.5)  U(x,t)  -  6(x-t)  +  R(xH)  ,   X  <  0  , 

and  U(x,t)  is  real  for  all  x  and  t. 

The  function  R(z)  has  the  property 

(3.^^)  R(z)  =  0  ,  z  <  0 

because  of  (2.1i)  and  the  regularity  of  r(k)  in  the  upper  half  plane.  The 
physical  meaning  of  (3.6)  in  connection  with  (3*5)  is  easily  understood, 
for  the  reflected  transient  should  not  appear  at  any  point  x  until  the 
appropriate  amoxrnt  of  time  has  elapsed  after  the  incident  transient  has 
struck  the  edge  of  the  reflecting  medium  at  x  ■  0. 

It  is  reasonable  to  expect  that  any  solution  U(x,t)  of  (3 '3)   can  be 
obtained  by  applying  a  linear  transformation  to  some  solution  U  (x,t)  of  the 
free  space  wave  equation 

2      2 

(3»7)  0 0  ,  Q 
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In  fact  it  would  not  be  unreasonable  to  suppose  that  there  exists  a  kernel 
function  K(x, z)  such  that 

00 

(3.8)  U(x,t)  -  U  (x,t)  +  I    K(x,z)  U  (z,t)dz  . 

o       1  o 

^-00 

Now,  if  U(x,t)  is  a  right  moving  transient  in  the  sense  that 

(3.9)  U(x,t)  -  0  ,  X  >  t  , 

then  it  would  seem  natural  to  suppose  that  U(x,t)  will  depend  only  on  U  (x,t) 
at  those  values  of  the  space  variable  x  vrtiich  the  non  zero  part  of  the 
traveling  wave  would  have  had  time  to  reach  at  the  time  t.  In  other  words, 
one  might  guess  that 

(3.10)  K(x,z)  =0,  z  >  X  . 
Assuming  (3.9)  and  (3,10),  (3.8)  becomes 

X 

(3.11)  U(x,t)  -  U  (x,t)  +  (   K(x,z)  D  (z,t)dz  . 

0         I  o 

^-00 

If  U  (x,t)  is  the  particular  solution  of   (3.7): 

(3.12)  U^(x,t)  -  5(x-t)   +  R(x+t)  , 
then   (3.11)  becomes 


(3.13) 


irtiere 


X 

U(x,t)  =  5(x-t)   +  R(x+t)   +  K(x,t)')^(x-t)   +   f     K(x,z)^(x-z)R(z+t)dz  , 

^(z)-/°'        ^<° 
^^  1 1  ,        z  »  0  . 


Assuming  that 

(3.1h)  K(x,t)  -  0  ,   t  <  -  X  , 

it  follows  from  (3.13)  that 


U(x,t)  -  U^(x,t)  ,   x<  0  , 
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It  can  be  shown  directly  that  there  exists  a  function  K(x, z)  satisfying 
(3.10)  and  (3.lU)  such  that  (3.11)  holds  l-^-l .  Applying  (3.m)  to  (3.11)  one 
obtains 


(3.15)  U(x,t)  »  < 


r  X 

^(x,t)  +    K(x,z)  U  (z,t)dz  ,  X  >  0 

0  J  0  ~" 

-X 


U  (x,t)  ,  X  <  0  . 


After  substituting  from  (3.15)  into  (3.3),  using  (3.7),  it  is  easy  to  verify 
that  (3.15)  provides  a  solution  U(x,t)  of  (3.3)  if  K(x,z)  satisfies 

2    2 

(3.16)  i|  .  1^  -  V(x)K  -  0  ,   z  <  X  , 

dx    3z 

subject  to  the  conditions 

(3.17)  K(x,-x)  =  0,  2^^^  =  V(x)  .* 

Since  it  is  known  from  the  theory  of  the  linear  second  oixier  hyperbolic  partial 
differential  equation  that  the  Gauchy  problem  defined  by  (3.16)  and  (3.1?) 
has  a  solution,  the  existence  of  K(x, z)  satisfying  the  required  conditions 
follows  and  thus  the  expression  (3.15)  for  U(x,t)'-^-J. 

Now  consider  (3.15)  with  U  (x,t)  given  by  (3.12).  Assume  that  the  cor- 
responding U(x,t)  is  a  right  moving  wave;  i.e., 

(3.18)  U(x,t)  =  0  ,   X  >  t  . 
Then  (3.15)  implies 

X 

(3.19)  0  -  R(x+t)  +  K(x,t)  +  I   K(x,z)R(z+t)dz  ,   x  >  t  . 

-00 


Conditions   (3.17)  and  continuity  of  K(x,z)  in  the  region  z  <  x  imply 

kf     V(x)dx:  ,       x  >  0 
K(x,x)  '  \jo 

0  ,       X  <  0 
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Since,  according  to  (3*6),  R(z+t)  is  zero  for  z<-t  the  lower  limit  in 
the  integral  of  (3.19)  can  be  replaced  by  -t,  and  it  follows  from  (3»19) 
that 

(3.20)  K(x,t)  -  0  ,   t  <  -X  . 

Equation  (3.19)  is  a  Fredholm  integral  equation  for  K(x,t)  considered 
as  a  function  of  t  (with  x  regarded  as  a  parameter)  with  the  kernel  R(z+t) 
over  the  range  -co  <  t  <  x.  When  R(z+t)  is  given,  e.g.,  by  experimental  data 
through  (3.h),  the  integral  equation  (3.19)  can  be  solved  for  K(x,t)  at  every 
value  of  x.  The  function  V(x)  can  then  be  obtained  from  the  second  relation 
in  (3.17). 

The  existence  of  a  solution  K(x,t)  of  (3.19)  can  be  established  by  means 
of  a  standard  argument  in  the  Fredholm  theory:  one  simply  proves  the  uniqueness 
of  the  solution.   At  the  same  time  the  uniqueness  of  the  solution  K(x,t)  of 
(3.19)  can  be  used  to  establish  the  fact  that  K(x,t)  satisfies  (3.16)  when 
V(x)  is  defined  by  (3.17).  The  first  relation  of  (3.17)  follows  at  once  fran 
(3.19)  and  (3.6).  The  fact  that  K(x,t)  satisfies  (3.16)  follows  without  much 
difficulty  on  applying  the  differential  operator 

2     2 

—5-  -  — 5-  -  V(x) 
dx    dt'^ 


to  (3.19).  The  result  is 

f 

(3.a )  0  - 


2     2 

— 5" 5"  -  V(x)K 

dx    dt"^ 


X 

-00 


'^^-^-  V(x)K 

lax   at 


,R(z+t)dz,  X  >  t  . 


Equation  (3.21)  is  the  homogeneous  integral  equation  corresponding  to  (3.19) 

for  the  function 

2     2 
^  K   a  K    u/^'\v 
— 5" 5  -  V(x;K  . 

ax'^  at*^ 


See  Section  k   and  [5^  • 
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If  the  solution  of  (3.19)  is  unique,  the  solution  of  (3.21)  must  vanish: 

^.L^    .     V(x)K  =  0 
dx  dt 

as  required  by  (3.16). 

h.       The  Case  in  Which  r(k)  is  a  Rational  Function 

If  the  reflection  coefficient  r(k)  is  known  for  all  real  values  of  k,  or 
even  for  just  the  positive  real  values  because  of  (2.5),  the  function  V(x) 
and  any  solution  u(x, k)  of  (1.3)  can  be  determined  as  indicated  in  Section  3. 
In  principle,  it  is  then  possible  to  find  the  other  elements  of  the  scattering 
matrix  in  terms  of  r(k).  As  a  practical  matter,  however,  the  relationship 
thus  exhibited  between  the  elements  of  the  scattering  matrix  is  quite  com- 
plicated and  indirect.  On  the  other  hand,  if  the  given  r(k)  is  a  rational 
function  of  k,  an  explicit  relation  between  r(k)  and  t(k)  or  between  r(k) 
and  p(k)  can  be  found  as  well  as  a  closed  form  expression  for  V(x).  These 
relations  between  the  scattering  matrix  elements  will  now  be  discussed. 

Consider  the  general  form  which  r(k)  can  take  if  it  is  a  rational  function 
of  k: 

(h.l)  r(k) ^ , 

TT  (k-X  ) 
i-1    ^ 

where  r  ,  the  u.  and  the  X.  are  complex  constants,  and  the  integers  ra  and  n 
0      1         i 

satisfy  the  inequality 
(U.2)  m  <  n  +  2  . 

The  condition  (U.2)  is  necessary  if  a  piecewise  continuous  V(x)  is  desired, 
i.e.,  if  r(k)  is  not  o(i)  for  large  k,  V(x)  will  contain  a  part  which  is 
improper,  e.g.,  a  Dirac  delta  function.  This  can  be  seen  by  examining  the 
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method  of  Section  3.  The  complex  constants  r  ,  X.,  [i.   must  be  such  that  the 
reality  condition  on  r(k:)  given  by  (2.5)  is  satisfied.  Also,  to  comply  with 
condition  (3.1)  which  depends  on  (3.6)  and  to  prevent  the  occurrence  of  proper 
eigenfunctions'-  -J  r(k)  must  be  regular  in  the  upper  half  of  the  complex  k 
plane.  In  other  words,  the  constants  X,  are  all  in  the  lower  half  plane: 

(Li,3)  Im  X  <  0,  i  -  1,2,  ,,.,  n  , 

Finally,  it  is  convenient  to  impose  the  condition 

(U.U)  |r(k)I  <  1,  for  all  real  k  . 

Condition  (b.'-i)  is  slightly  stronger  than  the  necessary  one 

(U.5)  |r(k)|  <  1,  for  all  real  k, 

implied  by  (2,8),  but  it  is  not  essentially  more  restrictive  since  the  r(k) 
satisfying  (U.5)  can  be  multiplied  by  a  parameter  e  such  that 

0  <  e  <  1 
and  the  correct  result  obtained  as  a  limiting  case  when  e  approaches  one.  In 
the  same  spirit,  it  will  be  assimed  that  the  numbers  x,  ,  which  will  be  intro- 
duced later  as  zeros  of  l-r(k)r(-k)  in  the  upper  naif  plane,  are  all  distinct, 
the  case  of  multiple  zeros  being  treated  as  a  limiting  case  in  which  two  or 
more  zeros  coalesce. 

If  some  a  priori  assimptions  are  made  about  t(k),  by  a  standard  argument 
the  connection  between  t(k)  and  r(k)  can  be  given  even  for  more  general 
reflection  coefficients  r(k).  Suppose,  in  fact,  it  is  assumed  that  t(k)  has 
no  zeros  or  singularities  in  the  upper  half  of  the  complex  k  plane.  For  real  k, 
by  (2.8)  and  (2.5): 


(U.6)  !t(k)I  -Vi-|r(k)|2  =  V^-r(k)r(-k)  . 
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The  function 

(h.7)  log  t(k)  =  log|t(k)|  +  i  e(k)  , 

k  real,  is  analytic  on  the  real  axis  and  in  the  upper  half  plane  according  to 
the  assumptions  made  about  t(k).  The  phase  e(k)  of  t(k)  is  an  odd  function 
of  k  and  the  modulus  |t(k)|  is  an  even  function  of  k  when  k  is  real  according 
to  (2.5).  Also  by  (2.U) 

(U.8)  lira    log  t(k)  =  0  , 

k  — >  CO 

if  the  limit  is  approached  on  the  real  axis  or  in  the  upper  half  plane.  By 

Cauchy's  integral  theorem,  integrating  over  an  infinite  semi-circle  in  the 

upper  half  plane  and  over  the  real  axis,  and  by  (U.8) 

(1.9)  lo,  tW  .  ^    f   IS^^  a^ , 

-00 

k  in  the  upper  half  plane.  Moreover,  since  Ic  is  outside  the  closed  contour 

-CD 

k   in  the  upper  half  plane.  On  taking  the  complex  conjugate  of  both  sides, 
(Li.lO)  becomes 

(U.ii)       o.-^fl^^.^, 

-00 

k  in  the  upper  half  plane.  Subtracting  (U.ll)  from  ()4.9)  gives  the  result: 

-00 

k  in  the  upper  half  plane,  where  (U.?)  and  (U.6)  have  been  used.  Equation  (U.12) 
provides  the  general  relation  between  t(k)  and  r(k)  when  t(k)  is  regular  and 
has  no  zeros  in  the  upper  half  plane.  If  these  conditions  are  allowed  to  be 
violated,  the  zeros  and  singularities  of  t(k)  must  be  known  in  order  to  specify 
a  relationship  between  r(k)  and  t(k)  by  means  of  an  argument  such  as  the  one 
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just  used,  based  entirely  on  the  general  properties  of  the  scattering  matrix 
elements  given  in  Section  2.  Apparently,  if  this  information  were  lacking, 
it  would  be  necessary  in  some  way  to  include  more  of  the  relationship  of  the 
scattering  matrix  to  the  differential  equation  (1.3 )• 

If  (U.12)  is  applied  to  the  r(k)  given  by  (b.l),  an  explicit  formula  for 
t(k)  can  be  obtained.  The  evaluation  of  the  integral  in  (U.12)  can  be  carried 
out  by  considering  first 

00 
-CD 

k  in  the  upper  half  plane.     After  an  integration  by  parts    (Li.13)  becomes 

d  1   r  a5'[^^^)^^-^] 

(H.IU)  dkl°g*(^)  =  -2^J      (o^-k)^l-r(o^)r(-a')}  ^"^  * 

-CD 

k  in  the  upper  half  plane.  The  integrand  of  the  integral  on  the  right  of  (h.lU) 
is  now  a  rational  function  of  a  and  can,  therefore,  be  evaluated  by  the  method 
of  residues.  The  result  is  a  rational  function  of  k  equal  to  the  derivative 
of  log  t(k).  The  integration  of  that  function  is  easily  carried  out  and  the 
constant  of  integration  determined  by  the  condition  (h.B).  The  result  is 

TT  (k*x^) 

(U.15)  ^(^^"ir ' 

TT  (k-x  ) 

i=l    ^ 
where  the  n  complex  constants  x.  are  the  n  roots  in  the  upper  half  plane  of  the 
equation 

(U.16)         l-r(k)r(-k)  «  0  . 

That  there  are  exactly  n  roots  of  (U.l6)  in  the  upper  half  plane  is  clear  from 
the  invariance  of  (U.16)  when  k  is  replaced  by  -k  and  from  condition  (U.U). 
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The  numbers  x.  are  assumed  to  be  distinct,  i.e.,  it  is  assumed  that  there  are 
1 

no  multiple  roots  of  the  equation  (U.l6).  If  this  is  not  so,  r(k)  can  be 
modified  slightly  so  that  the  multiple  roots  of  (I-1.I6)  split  up  into  simple 
roots,  and  one  can  consider  the  original  case  as  a  limiting  form  in  which  two 
or  more  roots  of  (U.16)  coalesce.  It  is  obvious  that  all  conditions  on  t(k) 
are  satisfied  by  (U.l^)  in  any  case,  even  if  one  ignores  its  derivation. 

The  phase  law  of  (2.8)  provides  the  means  for  obtaining  an  expression  for 
the  function  p(k)  once  t(k)  is  assumed  to  be  given  by  (Li,l5); 

.    m        n 

(U.I7)  p(k)  ' ^^^:-- ^ . 


TT  (k-\)  TT  (k-x.) 
i=l    ^  j=l    "^ 


Notice  that  p(k),  unlike  r(k),  has  poles  in  the  upper  half  plane:  the  x..  Since 
no  proper  eigenfunctions  can  exist  it  is  clear  from  Section  3  that  the  function 
V(x)  does  not  vanish  identically  for  all  values  of  x  larger  than  some  finite 
number;  i.e.,  V(x)  must  be  different  from  zero  in  the  positive  x  direction  at 
arbitrarily  large  values  of  x. 

Unfortunately  the  results  of  (U.15)  and  (U.l?)  cannot  be  trusted  on  the 
basis  of  the  argument  given  for  their  derivation  because  of  the  a  priori 
assumptions  made  on  the  t(k)  that  it  has  no  zeros  or  singularities  in  the 
upper  half  plane.  Instead  of  attempting  to  show  somehow  that  these  assumptions 
are  Justified  by  properties  of  the  differential  equation  (1.3)  there  will  be 
given  an  independent  derivation  which  will  at  the  same  time  provide  a  closed 
form  expression  for  V(x). 

The  relations  ()j,  15)  and(>i.l7)  are  demonstrated  by  starting  with  p(k)  as 
given  in  (^.I7)and  treating  it  as  the  reflection  coefficient  on  the  left  for 
the  differential  equation  (1.3)  with  x  replaced  by  -x.  The  reflection  coefficient 
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on  the  right  in  this  case  will  be  r(k).  Because  V(-x)  extends  to  infinity  on 
the  left  the  solution  u(-x,k),  which  can  be  obtained  from  p(k:)  by  the  method 
of  Section  3,  will  be  found  to  have  a  comparatively  simple  form.  If,  instead, 
r(k)  were  used  to  obtain  the  solution  u(x, k)  corresponding  to  V(x),  which 
vanishes  for  negative  values  of  x,  it  would  be  found  that  u(x,k)  is  much  more 
complicated  in  form.  The  demonstration  can  be  carried  out  by  showing  that 
u(-x,k)  and  its  derivative  at  x  =  0  match  the  assumed  form  of  u(-x,k)  and  its 
derivative  given  in  terms  of  r(k)  at  positive  values  of  x« 
Temporarily  it  is  convenient  to  define: 

^,    m        n 
(-l)'"-"*^  r^  IT  (k+^i)  TT  (lc*Xj) 

(U.18)         ^(k)  - ^^L_ tl 


TT  (k-x.)  TT  (k-x  ) 
i=l      ^   j=l      -^ 


and 


(U.19)         t(k)  -  ^=^ 


IT  (k-x.) 


TT  (k-X^) 

3=1 


Lemma:  for  all  real  k 


$(-k)  =  ?(k),  £(-k)  "  t(k)  , 
(U.20)         Ir(k)!  =  |$(k)|,  |r(k)|2  *   |t(k)|2  -  1  , 


F(k)  t(k)  +^(k)^(k)  -  0  . 
Proof:  Consider 


(U.21)         P,(k)  -  ?(k)^(-k)  fr(k-X.)  TT(-k-X,) 
^  i-1    ^  J^    '^ 
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and 


(U.22)  P,(k)  -  {l-r(k)r(-k)]  TT  (k-X,  )  JT  ("k-Xj  . 

^      ^  i=l    ^  j=l     J 

From  (U.l)  and  (U.19)  it  can  be  seen  that  P- (k)  and  Pp(k)  are  each  polynomials 
in  k  of  order  2n.  The  roots  of  P, (k)  are  the  2n  numbers  +  x.,  i=l,  ...,  n  . 
By  the  definition  of  the  x,  in  connection  with  (U.16)  and  the  fact  that 
replacing  k  by  -k  in  ()-i.l£)  leaves  the  equation  unchanged  it  is  obvious  that 

the  numbers  +  x.  ,  i=l,  ...,  n  are  also  the  roots  of  Pp(k).  After  dividing 

n        n 
(U.21)  and  (U.22)  by  JJ   (k-X  )  7T  ("1^~X . )  and  allowing  k  to  approach  oo  it  is 

i=l    ^  j=l     ^ 

found  that  each  of  the  expressions  approaches  one.  The  conclusion  is 

(U.23)         t(k)t(-k)  =  1  -  r(k)r(-k)  . 

It  follows  from  ()4.23),  (U.U)  and  (2.5)  that  t(k)t(-k)  is  real  and  positive  for  real 


By  (h.l)  and  (2.5),  for  real  k, 

F*fr  (k-il^) 


(U.2b)  r(-k)  -    ^'^ 


fr  (k-i.) 

i-l    ^ 


Analytic  continuation  implies  the  validity  of  (U.2U)  for  all  k.  From  ()4.l), 

(U.23)  and  (U.2U)  it  can  be  seen  that  the  2n  poles  of  ^(k)'?(-k)  are 

X . ,  i"l,  .,.,  n  and  X. ,   1=1,  ,..,  n.  Each  X.  is  distinct  from  all  the  T. 

since  the  X.  are  given  in  the  lower  half  plane.  It  can  then  be  seen  that  the 

2n  roots  of  (U.16)  consist  of  the  2n  numbers,  +  x.,  i-1,  ...,  n  as  well  as 

the  2n  numbers  +~xT,  i»l,  ...,  n.  Of  course  (U.l6)  has  at  most2n  distinct 

roots  altogether;  by  the  assumption  of  this  article,  in  fact,  it  has  exactly 

2n  distinct  roots.  By  (U.23),  therefore,  t(k)^(-k)  has  as  zeros  the 

±  X.,  1=1,  ...,  n  as  well  as  the  ±  x. ,  i"l,  ..*,  n.  Thus,  the  rational  function 
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n        n 

tt  (k+x  )  TT  (k+x  ) 

t(k)t(k)  =  ±^ 1^ 

TT  (k-x.)  T[  (k-x  ) 
i=i      ^   j=i      ^ 

A     A 

has  the  sajne  poles  and  zeros  as  t(lc;t(-k),  and  both  expressions  approach  one 
as  k  approaches  oo  .   It  follows  that 

(U.25)  ^(-k)  =  flkl,   k  real. 


From   (U,18)  and   ()4,19): 

(U.26)  p(i<)   -?737r~ 


The  relations  (li.20)  of  the  lemma  follow  at  once  from  ()4.25),  (Li.26)  and 
(l;.23). 

The  argtiments  just  presented  concerning  the  zeros  of  t(k)t(-k)  demonstrate 
at  the  same  time  that  the  n  zeros  of  t(-k)  are  the  x.,  i=l,  ...,  n  and  the 
-X.,  i=l,  ...,  n.  It  follows  from  (U.26)  that  the  n  poles  of  p(k)  in  the 
upper  half  plane  are  the  x.,  i=l,  ...,  n  and  the  -x.,  i=l,  ...,  n,  since  the 
X.  are  all  given  in  the  lower  half  plane.  By  the  Mittag-Leffler  theorem  the 
function  p(k)  can  be  written  in  the  form 


(U.27) 


n   Pv    JL  Pn+v 


V=l     V    V=l     V 


By  considering  the  analytic  continuation  of  0(k),  k  real,  and  comparing  with 
p(-k),  using  (U.20),  it  is  observed  that  (U.27)  can  be  expressed  as  the  sum  of 
pairs  of  terms  of  the  form 

Pi  P4  -Pi  pT 

T— =-  +  T-^  =     r— = T— ^  .   k  real  . 

k-oT       k+oT         Ic+oT       Ic^oT  >   "•  ^=  -^   • 


i  i  i       "^  ^i 

It  follows  at  once  then  that 


(U.28)  pj     -     -p^ 
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5.   V(x)  When  r(k)  is  a  Rational  Function 

Consider  the  problem  of  determining  V(x)  in  (1.3)  when  the  function  Zi^) 
given  by  (U.18)  is  the  reflection  coefficient  at  x  =  0.   The  kernel  of  the 
integral  equation  (3.19)  is 

CD 

-ix  (z+t)     .   „ 
e   v^     ,  z+t  <  0  , 


(5.1)  R(x+t)  -  ^  f  ^(k)e"^^^^** W  =  i  r  P^^ 


-00 

where  the  p  are  the  residues  of  p(k)  at  its  poles  in  the  upper  half  plane. 
From  the  lemma  of  Section  h   it  can  be  concluded  that  R(z+t)  is  real  for  real 
values  of  z  and  t. 

To  solve  the  integral  equation  (3.19)  it  is  convenient  to  assume  that 
when  X  is  negative  K(x,t)  has  the  form 

(5.2)  K(x,t)  "  Tf  (x)e"%"^  . 

With  this  assumption  the  integral  equation  (3.19)  takes  the  form 

n   -i(x  +x  )x 

(5.3)  -p  r  ^ JL_Ji_  f  (x)  +  f  (x)  +  ip  e"^V  =  0,  V  =  1,  ...,  n, 

1J.=1     V   p.      ^ 

which  holds  when  x  <  0. 

The  system  (5.3)  has  a  unique  solution  for  the  f  (x).  This  fact  can  be 
proved  in  exactly  the  same  way  that  the  uniqueness  of  the  solution  of  the 
integral  equation  (3.19)  can  be  proved^  -I . 


It  will  be  assumed  throughout  that  the  poles  of  ^(k)  in  the  upper  half  plane, 
the  X.,  are  distinct.  If  there  are  multiple  poles,  the  reflection  coefficient 
is  to  c)e  replaced  by  one  for  which  the  x.  are  distinct,  the  poles  which  were 
formerly  multiple  being  replaced  by  simple  poles  very  near  each  other.  The 
case  in  question  then  becomes  a  limiting  case  in  which  certain  simple  poles 
coalesce  to  form  multiple  poles. 
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Multiply  the  homogeneous  system  of  algebraic  equations  corresponding  to 
(^•3)  by  e  ^  and  sum  over  v: 


(5.h)     Z.  - 


n  o  e"^(^;■^^x)^   . 


-     X  +x 
V,p,»l      V   |J. 


^     v=l 


<  X  <  0 


Equation  (5»U)  can  be  written 


'  dz   E  ip  e-i^'(**^^  ^  e-^V  f°(x)l  +  X:  e"^""^*  f° 

U    r-1  ^^"1        ^  J   v-1 


(x) 


(^.5) 


X    r'   OD 

dz  , 


-CD 


^  f  $(lc)e-^'^(^-^*)dk  J:  e-^V  f°(x)  +  f  e-i'^vt  ^c 


(x)  -  0  , 


for  x+t<0,t<x.     On  setting 


(5.6) 


n 


-ix»,t  -o. 


U"(x,t)  -  Y:  e'^^^^  f°(x)  , 
v-1  ^ 


equation  ($*$)   becomes 

X 

(5.7)  I  R(t+z)U°(x,z)dz  +  U°(x,t)  -0,   x  +  t<0,  t<x. 

-CD 


Equation  (5.7 )  can  be  written 


(5.6) 


lim 
N  — >  CD 


-OD 


N 
1  (    A/,x  -ik(z+t) 

-N 


N         > 
dk  *  ^i  I  e-^^^^-*W  U°(x,z)d 


z  »  0, 


x  +  t<0,  t<x, 
since  U  (x,t)  is  obviously  absolutely  integrable  in  the  interval  -oo  to  x  . 
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Using  the  lemma  of  Section  k   it  can  be  verified  in  a  straightforward 
way  that 


($.9) 


N  f 
-N  L 


-ik(zH)    ^  ^-ik(z-t)(^^j^ 


|l.|^(k)|2[eik(t-z)^^  .  f  jeil^  .  $(k)e-^^  il^'il^  .  ^Me'^'^U   . 


-N 


After  substituting  (5.9)  into  (5.B),  multiplying  by  U  (x,t)  and  integrating 
over  t  from  -co  to  x,  there  results: 


lim 
N  — >  00 


,l-|p(k)|' 


r  X 


>  (       X 


ttO  /  .  \  -ikt ,. 
U  (x,t)e    dt 


-00 


,,o,   X  -ikz, 
U  (x,z)e    dz 


-00 


,dk 


(5.10) 


+   lim 
N  — >  00 


-N 


r  X 


-00 


e   +p(k)e 


U''(x,t)dtl 


-00 


ikz  A  /,  \  -ikzl;7T 


e   +p(k)e 


U  (x,z)dzLdk  «  0  , 


X  <  0  . 


According  to  (U.20),  the  lemma  of  Section  h,   conservation  of  energy  for 
p(k)  and  t(k)  holds  when  k  is  real;  i.e.,  |p |  <  1  when  k  is  real.  It  follows 
that  each  term  on  the  left  of  (5.10)  is  positive  for  any  function 


X 

/. 

-00 


U  (x,t)e    dt  unless  it  is  identically  zero  when  x  <  0,  and  so  the 


equation  (5.10)  implies  that  U°(x,t)  -Oifx+t<0.  From  (5.6),  then 


(5.11) 


f°(x)  =  0  ,   x  <  0  . 


Thus,  the  only  solution  of  the  homogeneous  system  (5.U)  is  the  trivial  solution 


f  2  0,  and,  therefore,  the  system  (5.3)  has  a  unique  solution. 
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For  a  reflection  coefficient  r(k:)  given  by  (U.l)  satisfying  the  general 
conditions  of  Section  2  as  well  as  the  special  conditions  (U»2),  (U.3)  and 
(U.'i)  relations  for  the  other  elements  t(k)  and  p(k)  of  the  scattering  matrix 
can  be  given.  In  fact, 

Theorem  1;  The  functions  p(k)  given  by  (U.lB)  and  t(k)  given  by  (b.l9) 
are  the  scattering  matrix  elements  p(k)  and  t(k)  associated  with  r(k)  given  by 

()4.1). 

Proof;  Suppose  -x  is  substituted  for  x  in  the  differential  equation 
(1.3),  30  that  V(-x)  replaces  V(x)  and  u(-x,k)  replaces  u(x,lc).  Set 

(5.12)  W(x)  =  V(-x),  v(x,k)  =  u(-x,k)  . 

The  reflection  coefficient  on  the  left,  r(k),  corresponding  to  V(x)  becomes  the 
reflection  coefficient  on  the  right  corresponding  to  W(x).  The  reflection 
coefficient  on  the  left  corresponding  to  W(x)  is  p(k),  and  the  transmission 
coefficient  on  either  side  is  t(k). 

Let  v(x,k)  be  the  solution  which  on  the  right  is  given  by 

(5.13)  v(x,k)  -  e"^'^  +  r(k)e^^  ,  x  >  0  , 
and  on  the  left  by 

(5.m)         lijn   ]v(x,k)  -  t(k)e"^'°'[  =  0  . 

X  >  -00  L  ) 

Let  v(x,k)  be  the  solution  which  is  determined  on  the  left  by 

(5.15)  llm    |v(x,k)  -  t(k)e"^^r  =  0  . 

X  — >  -00  [  J 

To  prove  the  theorem  it  is  only  necessary  to  prove  that 

(5.16)  v(x,k)  =  e"^""  +  r(k)e^^  ,   x  >  0  , 


-2U- 
or  equivalently  that 

(?.17)         v(o,k)  =  l+r(k)  and  gl  ^  =  ik  <   r(k)-l  > 


{ 


From  (5.15),  by  taking  the  Fourier  transform  of  (3.11)  and  substituting 
(5.2),  the  following  expression  for  v(x,k)  is  obtained: 


v(x,k)  =  t(k)e     +  it(k)e     2_  ^  ^u" 

v=l   V 


,   X  <  0 


The  relations  (5»17)  can  now  be  written 


(5.19) 


t(k)  +  it(k) 


n  f,.(0) 


V=l   V 


(k) 


and 


(5.20) 


A    A    n  A     n   f'(0)      , 

-ikt(k)+t(k)  21  ^^(0)   +  it(k)  XI  -^  •  ik  {  r(k)-l 
v»l  V^   V  "^ 


In  order  to  demonstrate  (5.19)  it  is  convenient  to  define  the  complex 


numbers  g  by  the  equation 


(5.21) 


f  ^ 

V"l    V 


-i 


r(k)  ^   1  _ 


^(k)   ?(k) 


The  expression  on  the  left  of  (5.21)  is  a  legitimate  representation  of  the 
rational  function  of  k  on  the  right  by  the  Mittag-Leff ler  theorem,  and  the 
coefficients  g  are  given  by 


(5.22)  g  -  -i  lim   (k+x  )  . 
k  ->-x    ^ 

V 


it-ic) .  JL_ .  u  .  -1 


'    V   .  1 


t(k)   ^(k) 


^^ 


,   v»l,  ,,.,  n 
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-V 


where  the  residue  of  jvm  i  at  the  pole-x   has  been  designated  by  ^     , 
(5.3) 

(5.23) 


From 


^(0)=p}t    JL_-il  v  =  l,  ...,  n. 
[^  M=l  V  tx    J 


The  numbers  g  of  (5.22)  will  satisfy  the  system  of  equations  (5.23)  when 
substituted  for  the  numbers  f  (O).  To  see  this,  consider 


(5.2U) 


n   g. 


-  ii  =  -ip 


f^(\) 


M^(x^)    ^(%) 


,  ,  V  =  1,  ...,  n  , 


by  (5.21).  The  right  side  of  (5.2U)  can  be  written 


lim    -  i(k-x  )p(k) 


r(k)    1  , 


(5.25) 


k  — >  X      ^  lt(-l<:)        t(-k) 


-1. 


'    V   .   1 


-V     -V 


■^ —   "^  ^\>   ^  =  1»  .•»>  n  , 


making  use  of  the  lemma  relations  (I4.20).  From  (5.22),  (5.2U)  and  (5.25): 


(5.26) 


g  "  0 


n   g 

^  X  +x 

H.=l  V  |j. 


V  »  1,  ,,,,  n; 


i.e.,  the  g^  satisfy  the  sjrstem  of  equations  (5.23). 

Since  the  system  (5.23)  must  have  a  unique  solution 


(5.27) 


^v^°^  "  ^v  '   ^  "  ^*  •••»  "  • 


After  substituting  from  (5.2?)  into  (5.21)  the  relation  (5.19)  follows, 
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To  demonstrate  the  remaining  relation  (5«20)  it  is  necessary  to  consider 
the  derivative  of  (5.3)  at  x  =  0: 

n   f'(0) 
(5.28)    X  n  +  f  ■  .(0)  +  ip,,  y  f  ,(0)-D  y    -ii—  =  0,  V  =  1,  ...,  n  . 


M- 


>.28)      x^p^  .  f  ^(0)  .  ip^   Z  f,(o)-Pv  ^,  rW 

ti=l  |J.=1   V   I 

Let  the  numbers  h  be  defined  by  the  equation 

(5.29)    it(k)  f^  ^  =  ik  {r(k)-l)  +  ikt(k)  -  t(k)  ^  f  (0) 
V=l    V  v=l 

From  (5.21),  (5.2?)  and  (U.2)  it  follows  that 


n 

(5.30)         Z  ^v^°^  '  ^^      ^^  ' 
v»l        k  >  CD 


A 


t(k)-l 


Relation  (5.30)  is  consistent  with  (5.29)  as  can  be  seen  by  allowing  k  to 
approach  oo  in  (5.29)  Since  the  poles  on  the  left  side  of  (5.29)  are 
the  same  as  those  on  the  right,  (5.29)  is  entirely  consistent.  The  numbers 
h  are  then  given  by 

(5.31)   h^=  li.  (»c^x^)^%^;^^^^-^  -  i  t  ^M  r-  ^'7'^"''>-  °  '"'•> "  • 

Also,  from  (5.29): 

(5.32)  ^  \L.  =  ^  ^—^ }L_J  +  i  T-  f  (0),  V  "   1,  ...,  n  . 

|i»l   V  [i.  ''  v'  ^"^ 

After  substituting  from  (5.32)  and  (5.31)  into  the  expression 

(5.33)  p,f  ^r^      h^ 

p,=l    V   |i 
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there  results 


Pv\^K^\)-l}   .    n         X  {r(-x  )-l] 

|i=l  ^ 


T(x  ) 


'wy 


(5.3U) 


lim    (k-x)p(k)kllii^H^isll2: 


k  — >  X 


lira   (k-x  )k 
k  — >  X    ^ 

V 

n 


T(k7 


''  H=l   ^^       k  -»  X      ^   ^HO 


^ 


-^'li^-'^'^im! 


*^^^  l.'^'°\  i^«;""v'M|gi  -  ^,  ■  «^„  ^  ip„  Z^  r^(o),  „  .  ,, 


.,  n. 


using  the  relations  (U.20)  of  the  lemma.  The  identity  of  the  expressions  (5.33) 
and  (5.3U)  means  that  the  numbers  h  when  substituted  for  the  f *(0)  satisfy 
the  system  (5.28).  Since  (5.28)  has  a  unique  solution  it  follows  that 


(5.35) 


f^(0)  -  h^  ,  V  =  1,  ...,  n  . 


From  (5.35)  and  (5.29)  it  follows  that  the  relation  (5-.20)  is  satisfied. 
The  remarks  made  concerning  (5.19)  and  (5.20)  show  that 


(5.36) 


t(k)  »  t(k)  . 


The  fdct  that  the  scattering  matrix  is  unitary  implies  then  that 

(5.37)         p(k)  -  0(k)  , 

and  Theorem  1  is  proved. 

It  is  now  possible  to  obtain  a  closed  form  expression  for  the  function  V(x) 
corresponding  to  r(k).  This  result  is  contained  in 
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Theorem  2:  Given  a  reflection  coefficient  r(k)  of  the  form  (li.l)  subject  to 
the  general  conditions  on  the  scattering  matrix  elements  and  the  special 
condition  (h.2).  The  corresponding  function  V(x)  of (1.3)  is  then  given  by 

.2 

(5.38)  V(x)  =  -  -2^  logA(x),    x<0, 

dx 

where  A(x)  is  the  determinant  of  the  matrix  whose  elements  have  the  form 

_p^gi(xv+x^)x 

(5.39)  5   +3^-^        ,  V,  ^  =  1,  ...,  n  . 

**    V  M- 

Here  the  quantity  6  ^   is  the   Rronecker  delta 


(5.U0)  6, 


\IV 


0,  [i   ¥  V 

1,  ti  =  V 


and  the  x  are  roots  in  the  upper  half  plane  of  the  equation 

(5.U1)         1  -  r(k)r(-k)  -  0  . 

Proof:  If  the  solutions  f  (x)  of  the  system  of  equations  (5.3)  are  represented 
by  Cramer's  rule  as  ratios  of  determinants  and  are  then  substituted  into  (5.2) 
with  t  set  equal  to  x,  it  can  be  sdsn  easily  that 


(5.U2)         ^  e-i'V^  f  (x)  -  -  ll^gM^il  ,   x< 
(0,=!        ^ 


It  follows  at  once  from  (3.17)  that 

(5.143)        w(x)  -  -/^°gA(-x)  ^   ^  ^  0  ^ 

dx"^ 

and  from  the  definition  (5.12)  of  W(x)  the  theorem  follows, 
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Finally,  it  should  be  noted  that  two  linearly  independent  solutions  of 
the  differential  equation  (1.3)  are  given  by 

.-w  n   f  (-x)e^^^^ 

i^M  u(x,±k)  =  e-^'^  .  iei^'^  Y:    ^hrrv 

v=l    V^ 


6.  Conclusion 

For  r(k)  a  reflection  coefficient  at  x  =  0  which  is  a  rational  function  of  k, 
the  transmission  coefficient,  t(k),  exists,  is  regular  on  the  real  axis  and  in 
the  upper  half  plane  and  has  no  zeros  on  the  real  axis  or  in  the  upper  half 
plane.  The  expressions  for  the  scattering  matrix  elements  t(k)  and  p(k)  in 
terms  of  r(k)  are  given  in  Section  k   in  (U.19)  and  (U.18)  when  r(k)  is  given 
by  (U.l).  The  expression  for  the  corresponding  V(x)  is  given  in  Section  S   in 
(5»38).  The  closing  paragraphs  of  Section  5  following  the  proof  of  the  theorem 
provide  the  necessary  description  for  understanding  (5«38).  Two  linearly 
independent  solutions  of  the  original  differential  equation  (1,3)  are  given  in 

(5.UU). 
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